Scaling Solution For Small Cosmic String Loops 
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The equation governing the time evolution of the number density of loops in a cosmic string 
network is a detailed balance determined by energy conservation. We solve this equation with the 
inclusion of the gravitational radiation effect which causes the loops to shrink (and eventually decay) 
as time elapses. The solution approaches a scaling regime in which the total energy density in loops 
remains finite, converging both in the infrared and in the ultraviolet. 
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INTRODUCTION 

More than 30 years ago it was realized that a net- 
work of one-dimensional topological defects could form 
during a cosmological phase transition Such cosmic 
strings consequently evolve in a fashion determined by 
the expansion of the universe, intercommutation events 
and gravitational radiation. All these mechanisms act in 
such a way so as to yield a scaling regime, in which the 
network remains self-similar under rescaling, at least in 
its large scale properties. This has been convincingly es- 
tablished by now, both through numerical 0, d, Q, d, |g, 0] 
and analytic @, Q studies, and it is widely believed 
that the inter-string distance and the correlation length 
along a long string both scale. However, the existence 
or not of a scaling regime for the smaller scales (for ex- 
ample, the typical size of the small loops) still remains 
an open question (loj . It is this problem we wish to ad- 
dress in this short note. Understanding the properties 
of small cosmic loops would be a valuable achievement: 
the detectability of gravitational waves emitted by cos- 
mic strings 11 , T3, HI depends crucially on their knowl- 
edge, while the size of the smallest loops can also have 
some effect on microlensing [14| . To this end we will con- 
sider the number density of loops - it is the small scale 
structure that is responsible for their production and so 
the number density of small loops provides a measure of 
the former. Presently, we demonstrate that this quantity 
scales. 

There are many processes that can act to change the 
number density of loops with a given length I. Expansion 
of the universe causes the strings to stretch, but only on 
scales larger than the horizon size [H|]. Therefore loops 
essentially do not grow in time. 1 Self- intersections of 
cosmic strings can produce loops from long strings and 
also fragment loop s into smaller ones. This mechanism 
was studied in [16], taking as an input the simplified ana- 
lytic model developed therein, which we believe describes 
fairly accurately the evolution of the small scale structure 



1 Loops bigger than the horizon distance are regarded as long 
strings. 



on cosmic string networks. Of course, string intercommu- 
tation can lead to the absorption of loops back into the 
long string population but this process is strongly sup- 
pressed for small loops and so can be neglected. Finally, 
the coupling of matter in the form of cosmic strings to 
gravity means that the network radiates away part of 
its e nergy . This has been amply discussed in the litera- 
ture [UBS OS III 

and as a consequence loops shrink as 
time progresses and eventually disappear from the net- 
work. 

The equation governing the evolution of the number 
density of loops in a cosmic string network is a detailed 
balance determined by conservation of energy. In what 
follows we solve this equation, first ignoring gravitational 
radiation and then taking into account this effect. It is 
shown that the solution approaches a scaling regime and 
if we include GW emission the total energy density in 
loops remains finite, in particular converging in the ultra- 
violet (UV), even when the loop production function di- 
verges. Above the gravitational radiation scale the loop 
energy density decreases with the size of the loops as a 
power-law which is in good agreement with the numerical 
simulations of 0, 0] and implies the convergence of the 
total energy density in loops in the infra-red (IR). How- 
ever, below the gravitational radiation scale this quantity 
falls off with a lower power of the loop length in such a 
way that the energy density remains finite in the UV limit 
also. Finally, we comment on the curious fact that this 
sub-gravitational radiation regime seems to be apparent 
in the results of H, even though those simulations did 
not include gravitational radiation. 



PRELIMINARIES 

Several scales can be defined for a cosmic string net- 
work. Among them, the characteristic length of the net- 
work plays an important part. It is defined as the length 
scale L such that a typical volume L z of the network 
contains a length L of long strings. Another important 
concept is that of a scaling regime. We say that a quan- 
tity (with units of length) is scaling if it remains constant 
in units of the cosmological time t as the network evolves. 
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Thus, in a scaling regime we have L = "ft, for some con- 
stant 7, and the energy density in long strings becomes 
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where fx is the string tension. According to Ref. Q the 
constant 7~ 2 takes the value ~ 3 in a matter-dominated 
universe and ~ 9.5 in a radiation-dominated universe, 
whereas Ref. [6| provides the values 7~ 2 ~ 3 and 7~ 2 ~ 
11.5, respectively. 

Now, denote by 4rr(£,t) the number density of loops 
with size comprised between Z and I + dl. This quan- 
tity has units of (length) ~ 4 . Thus, if it ever reaches a 
scaling regime during the cosmological evolution, it must 
eventually approach the following form: 
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t- 4 f(i/t) 



(2) 



Borrowing notation from [f|, define the length of loops 
in units of the horizon size, a = l/d^. If the scale fac- 
tor takes the form a(t) oc t v , the horizon size can be 
expressed as dh = t/(l — v) and so the signature of a 
scaling regime in the loop number density is a solution 
of the form 
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da 
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(3) 



S (a) is known as the scaling function. 

In what follows we shall need expressions for the rate 
at which the loops shrink due to emission of gravitational 
radiation and for the rate of loop formation. The power 
radiated away by cosmic strings is typically given by P = 
TG/i 2 , where G is Newton's gravitational constant and T 
is a numerical constant of order 50 [H EI ED, [H HH . 
Since a loop of length I has energy /zZ, this implies that 
loops shrink at a rate 



dt = - TG » 



(4) 



A first attempt to analytically estimate the loop pro- 
duction from long strings [161 ] resulted in a divergent total 
energy density in small loops. There it was found that 
the average number of loops produced per unit time, per 
unit distance along the string and per unit loop length is 
given by 



d(M) _ c n 
dt da dl "Mt 
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(5) 



where c is the overall normalization of the loop produc- 
tion function. For a matter-dominated era \ = 0-25 and 
c = 0.042, while during the radiation epoch \ = 0.10 
and c = 0.121. These values of the normalization c are 
most likely to be over-estimates since fragmentation was 
not taken into account. Indeed, these values of c exceed 



the numerical results of [B| by factors of approximately 
15 and 40, while comparison with Q yields factors of 
3.5 and 85, respectively. The latter reference considered 
the loop production function and so their results relate 
more directly to the normalization c, whereas the former 
reference studies the loop number density. Note that 
even though these simulations might seem in conflict, 
the authors of Q emphasize that the range of lengths 
corresponding to the small loops is not yet scaling. Fur- 
thermore, the determination of the normalization from 
the scaling function is sensitive to the choice of the range 
where the fitting is performed in jBJ. So there is a possi- 
bility that there exists no real discrepancy between both 
simulations. 

For a consistent description of the network, equa- 
tion §5§ must be corrected at small scales. The cor- 
rection is provided by the smoothing due to emission of 
gravitational waves (GW) from cosmic strings and has 
been computed in [24j considering a discrete spectrum 
for the fluctuations. More recently an improved result 
performed in the continuum limit was obtained in [25j ]. 
where the structure along the strings was found to be 
smoothened on scales below 
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Hence, gravitational radiation introduces a natural cutoff 
for the divergent loop production function, even though 
it is believed to be unnecessary to achieve scaling [5j, [6( . 
However, we will disregard this effect in what follows and 
see how far we can get. 



EVOLUTION OF THE LOOP NUMBER DENSITY 
NEGLECTING GRAVITATIONAL RADIATION 

As discussed in [ll|, if gravitational radiation is ne- 
glected the number of loops within a comoving volume 
changes only due to loop production: 



d ( 3 dn 
dt\ a ~dl 
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Defining F(l,t) = ^jj(l,t) and inserting the power-law 
expression for the scale factor we obtain 



tF + 2,vF = 



1+2X 



7 2 Z 4 \t t 
which has the general solution 

'Z\ 2x " 3 , C (Z) 



F{l,t) 



7 2 (3j/- 1 - 2 X )t 4 \t 



t 3 " 



(8) 



(9) 



The function Co depends only on the variable Z. Not- 
ing that the inequality 1 + 2\ < 3f is satisfied both in 
the radiation- and matter-dominated eras, we conclude 
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that at late times (t 
proaches 



oo ) the loop number density ap- 



dn 



(l,t) 



7 2 (3^ - 1 - 2 X )t 4 \t 
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(10) 



Thus, expressing everything in terms of a and dh we in- 
deed find that the loop number density approaches a scal- 
ing regime, i.e. it takes the form ^ with 



S(a) 



cjl-v)- 1 - 2 * 
7 2 (3^-l-2x)' 



(11) 



This confirms that the loop number density does ap- 
proach a scaling regime without taking into account grav- 
itational radiation. The exponent 2 X — 2 = — p takes the 
values —1.5 in the matter era and —1.8 in the radiation 
era. This is in good agreement with the numerical results 



-0.08 
-0.07 



and p r 



1.60 



-0.21 



of [H who quote p ma t — 1-41" 
However, these simple power-laws become good fits to 
the data only above a physical length i c , which is iden- 
tified with the initial correlation length of the network. 
Furthermore, if the solution (fTTj) were valid over the full 
range of loop lengths the total energy density would di- 
verge in the UV since 



l n dl 



< 2 



S(a)da . 



(12) 



We will now show that including the process of shrinkage 
of the loops due to gravitational radiation changes the 
power-law below the gravitational radiation scale, thus 
yielding a convergent total energy density in loops. 



EVOLUTION OF THE LOOP NUMBER DENSITY 
INCLUDING GRAVITATIONAL RADIATION 

Inclusion of gravitational radiation into the evolution 
equation for the loop number density introduces an extra 
term on the left-hand side of equation (JSJ because loops 
shrink at a rate given by (UJ). Therefore, we now have 



t F + 2,vF - TG/iF' 
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Defining for convenience b = TG/i, the solution of the 
above differential equation may be written as 



F(l,t) 



:b{l + bt) 2 *- 4 ( bt 



V 



t 2 * 
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C{1 + bt) 

t 3u 
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where B represents the Euler incomplete beta function 
and C can be any function of the combination I + bt. 



By employing the Taylor expansion 

B(e;3v-l-2 X ,2 X -2) = /*" - + O(e) , (15) 

6v — 1 — 2 X 

we find that the first term in (fT"4l) behaves, for large I 
and fixed t, as ~ l 2x ~ 3 , whereas the second term goes 
like ~ C(l). Requiring that the energy density in loops 
converges in the limit I — > oo imposes that the general 
function C(x) decays faster than x~ 2 . Therefore, using 
the expansion 

fl(l-e;3i/-l-2x,2x-2) = 

7rcsc(27rx)r(3^ - 2x)e 2x ~ 2 

(31/ - 1 - 2 X )r(2 X - l)r(3 - 2 X )r(3^ - 1 - 2 X ) 

1 



o 



(16) 



the second term in ([T4|) is dominated by the first term as 
t — > oo: 



dn , 



cb(l + bt) 2 *- 4 



r 

x B 



f 2 x 



bt 



l + bt 



2 X -3v 



bt 



l + bt 



;3^-l-2 X ,2 X -2 .(17) 



Once again converting to the variables a and d^, we ob- 
tain a solution of the form ([3]) with 



S(a} = 
x B 



cb 2 *- 3 a fa + (l-v)b 



7 2 (1-^) 4 V (l-v)b 
(l-i/)6 



3v-4 



a+{l-v)b' 



3^-l-2 X ,2 X -2 
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Now we can use the series expansions (1151) and ((To)) to 
recover the limits for small and large loops. The separa- 
tion between these two regimes is set by the gravitational 
radiation scale TGfi, and we find for a 3> TGp the same 
result ITT)) , whereas for a <C TGfi 



S(a) 



7r(l-2 X )csc(2^ X )r(3^-2 X ) 
(31/ - 1 - 2 X )r(2 X )r(3 - 2 X )T(3v - 1 - 2 X ) 

(19) 



7 2 rG/^(l - i/) 2+2 x 



v 2 x-i 



Since < 2 X < 1 holds in both cosmological eras, the 
integral |(T2")) is manifestly convergent and so the total 
energy density in loops is finite, as desired. 



DISCUSSION 

We may now compare our results with those obtained 
numerically in [f|. We have already noted that the ex- 
ponent 2 X — 2 agrees well with the simulations for the 
larger loops. However, we have shown that taking into ac- 
count gravitational radiation has the effect of bending the 
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curves for a < TGfi so that the exponent then becomes 
2\ — 1. Indeed, Figure 3 of that reference does appear to 
show a certain range of the parameter a in which the scal- 
ing function behaves as such a power-law. At first sight 
this might seem intriguing since those simulations did not 
include the gravitational radiation process directly. How- 
ever, small loops behave like matter and so the expansion 
of the universe effectively shrinks the loops. Because the 
simulations only keep loops with sizes greater than a fixed 
fraction of the horizon, they are eventually removed from 
the game, hence emulating gravitational radiation 2 . The 
minimum counting size was a m i n — 1CP 5 so we should 
expect the sub-gravitational radiation regime to set in 
for comparable scales or smaller. A more accurate esti- 
mate, equating the lifetime of the loops determined by 
gravitational decay to the lifetime set by the minimum 
counting size, yields abend — j-r?;Ctmin for the scale at 
which the bending would occur. Nonetheless, it is cu- 
rious that this sub-gravitational radiation regime shows 
up also at early times in the simulations, when scaling is 
yet to be achieved. 

We also note that the inclusion of gravitational radi- 
ation in our equations leads to a finite energy density 
of loops, even though the loop production function used 
as an input diverges at small scales. This means that 
the rate at which loops are removed from the network 
is sufficiently high to balance the diverging loop forma- 
tion. Of course, this divergence must be eliminated in 
order to satisfy string length conservation. The gravita- 
tional smoothing lcngthscale (J6j> naturally introduces a 
UV cutoff on the loop production function and so this 
should be the typical size of cosmic loops, as advocated 
in 25j . It is comforting that this cutoff sets in much be- 
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fore the scale of validity of the effective theory, namely 
the thickness of the strings [26], is reached, as long as 
the value of the string tension is not too small. It re- 
mains to be seen if the the correct normalization of the 
loop production function can be reproduced analytically 
by methods similar to those used in 16]. This is left for 
future work. 
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